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ABSTRACT 


Motivated by recent THEMIS observations, this paper uses 2.5-D electromagnetic hybrid 
simulations to investigate the formation of Spontaneous Hot Flow Anomalies (SHFA) 
upstream of quasi-parallel bow shocks during steady solar wind conditions and in the 
absence of discontinuities. The results show the formation of a large number of structures 
along and upstream of the quasi-parallel bow shock. Their outer edges exhibit density and 
magnetic field enhancements, while their cores exhibit drops in density, magnetic field, 
solar wind velocity and enhancements in ion temperature. Using virtual spacecraft in the 
simulation, we show that the signatures of these structures in the time series data are very 
similar to those of SHFAs seen in THEMIS data and conclude that they correspond to 
SHFAs. Examination of the simulation data shows that SHFAs form as the result of 
foreshock cavitons interacting with the bow shock. Foreshock cavitons in turn form due 
to the nonlinear evolution of UFF waves generated by the interaction of the solar wind 
with the backstreaming ions. Because foreshock cavitons are an inherent part of the shock 
dissipation process, the formation of SHFAs is also an inherent part of the dissipation 
process leading to a highly non-uniform plasma in the quasi-parallel magnetosheath 
including large scale density and magnetic field cavities. 
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INTRODUCTION 


Collisionless dissipation processes at the bow shock result in reflection and/or 
leakage of ions into the upstream region forming the ion foreshock region (Asbridge et 
al., 1968; Greenstadt et al., 1968; 1980; Gosling et al., 1978; Paschmann et al., 1979; 
Bonifazi et al., 1980a,b). The ion foreshock is populated with a variety of ULF waves 
(e.g. Russell and Hoppe 1983; Le and Russell, 1992; Greenstadt et al., 1995) with wave 
vectors towards the sun but carried back by the solar wind in the opposite direction. Both 
observations and theoretical studies have also established the turbulent nature of the 
quasi-parallel shocks and the cyclic reformation of the shock front (e.g. Greenstadt et al., 
1977, 1993; Russell, 1988, Thomsen et al., 1988, Thomsen et al., 1990a, b; Burgess 1989; 
Thomas et al., 1990; Winske et al., 1990; Omidi et al, 1990; Scholer et al., 1993). This 
behavior is thought to be caused by the convection of upstream generated ULF waves 
into the shock. 


In an accompanying paper, Zhang et al. [2012] use THEMIS multi-spacecraft 
measurements to identify a new structure at the quasi-parallel bow shock named 
Spontaneous Hot Flow Anomaly (SHFA). SHFAs and Hot Flow Anomalies (HFAs) 
exhibit similar signatures in spacecraft time series data that consist of enhancements in 
density and magnetic field in the outer part and depletions in these parameters in the core 


which is also associated with increased temperature and deflected solar wind flow 


However, while HFAs form due to the interaction of solar wind discontinuities with the 


bow shock (e.g. Schwartz et al., 1988;1995;2000; Thomsen et al., 1986; 1988; 1993; 
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Paschmann et al., 1988; Thomas et al., 1991; Sibeck et al., 1998;1999;2000; Lin, 


1997;2002; Lucek et al., 2004; Omidi and Sibeck, 2007; Facsko et al., 2008; Eastwood et 


al., 2008; Jacobsen et al., 2009), SHFAs form in the absence of discontinuities. In the 
past, local and global hybrid (kinetic ions, fluid electrons) simulations have been used 
successfully to examine the formation and impacts of HFAs at the bow shock (e.g. 
Thomas et al., 1991; Fin, 1997; 2002 and Omidi and Sibeck, 2007). Motivated by SHFA 
observations, we have conducted an investigation of the quasi-parallel bow shock using 
global hybrid simulations. As we demonstrate here, simulations show the formation of 
copious structures at the quasi-parallel bow shock and foreshock whose time series 
signatures resemble those of SHFAs presented by Zhang et al. [2012]. The results 
indicate that SHFAs are an inherent part of the super-critical quasi-parallel shock 
dissipation processes and result in highly turbulent and non-uniform magnetosheath 
plasma. 


The structure of the paper is as follows. Section 2 describes the hybrid model used in 
this study while the simulation results are described in section 3. Section 4 provides a 
summary and conclusions. 


2. HYBRID SIMUFATION MODEF 

The main tool of investigation in this study is a 2.5-D (2-D in space and 3-D in 
currents and electromagnetic fields) global hybrid simulation model used extensively in 
the past (e.g. Omidi et al., 2004, 2005, 2006, 2009a, b; 2010; Omidi and Sibeck, 2007; 
Blanco-Cano et al., 2006a, b, 2009, 2011; Sibeck et al., 2008). In electromagnetic hybrid 
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codes, ions are treated as macro-particles and consist of one or more species (e.g., 
differing mass, charge, etc.) whereas electrons are treated as a massless, charge 
neutralizing fluid (see e.g. Winske and Omidi, 1993, 1996). 

The model consists of a dipole inside a sphere whose surface represents the 
ionospheric boundary. A solar wind type plasma with electron and ion betas (ratio of 
thermal to magnetic pressure) of 0.3 each and flow speed of 12 Va (Alfven speed) is 
uniformly loaded in the system except for the region inside the ionospheric boundary. 
This plasma is continuously injected from the left hand boundary throughout the whole 
run. The remaining boundaries remain open for the plasma to leave. Similarly, open 
boundary conditions are applied for the electromagnetic fields so that excited waves and 
turbulence in the system leave through these boundaries. The simulation box lies in the 
X-Z (noon-midnight meridian) plane with X along the solar wind flow direction (Sun- 
Earth line) and the magnetic dipole moment in the Z direction so that X corresponds to - 
Xgsm and Z corresponds to Zgsm. The simulation box extends 1500 ion skin depths c/co P 
(where c is the speed of light and co P is the ion plasma frequency) in the X and Z 
directions with cell size of 1 ion skin depth. The interplanetary magnetic field (IMF) lies 
in the X-Z plane and makes a cone angle of 10° with the X axis. To optimize the 
computational resources, the simulated magnetosphere is smaller (by a factor of ~5) than 
the Earth’s magnetosphere. On the other hand, the simulated plasma parameters and 
characteristic time and spatial scales such as gyroperiod, or ion skin depth are the same as 
in the solar wind and magnetosphere. This ensures that the simulations are capable of 
generating plasma and field values and characteristic scales that can be directly compared 


5 


141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 


to observations at the Earth’s bow shock. As demonstrated in our earlier studies, the 


physical processes occurring in smaller bow shocks and magnetospheres are similar to 
those at the Earth’s magnetosphere and much can be learned from these simulations 
including scaling properties of various magnetospheric processes (e.g. Omidi et al., 2004, 
2005, 2006, 2009a, b, 2010; Omidi and Sibeck, 2007; Blanco-Cano et al., 2006a, b, 2009, 
2011; Sibeck et al., 2008). 


3. FORMATION OF SHFAs 


Panel (a) in Figure 1 shows the plasma density (normalized to solar wind value) 
and magnetic field lines in a portion of the simulation domain. The quasi-perpendicular 
and parallel portions of the bow shock are labeled in this panel with the latter falling 
primarily in the southern hemisphere. Also labeled is the ion foreshock, upstream of the 
quasi-parallel shock, and the Foreshock Compressional Boundary (FCB) that separates a 
highly disturbed and turbulent ion foreshock plasma from a nearly pristine like solar wind 
that falls inside the ion foreshock (beam) boundary (see Sibeck et al., 2008; Omidi et al., 
2009b). Panel (b) in Figure 1 shows the density zoomed around the quasi-parallel shock 
and the ion foreshock. The latter includes regions of low density labeled foreshock 
cavitons. The presence of these structures was predicted by global hybrid simulations 


(Fin, 2003; Fin and Wang, 2005; Omidi, 2007) and confirmed in the ion foreshock 
(Blanco-Cano et al., 2009, 2011; Kajdic et al. 2010, 2011). Foreshock cavitons are about 
an Re (Earth radii) in size and are associated with drops in density and magnetic field in 
their core by as much as 50% or more and plasma and magnetic field enhancements in 
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their outer edge. They form as a result of the nonlinear evolution of ULF waves and are 
carried back by the solar wind towards the bow shock. As we show here, the interaction 
between foreshock cavitons and the bow shock is highly significant and an inherent part 
of the quasi-parallel shock dissipation processes. 


Although at any given time the structure of the quasi-parallel bow shock is highly 
turbulent, a closer examination reveals processes that occur at and upstream of the shock 
on a regular basis. An example of this is illustrated in Figures 2 and 3 that show the 
density and ion temperature (normalized to solar wind value) respectively at 4 different 


times (normalized to proton gyroperiod Q [ ) zoomed around the quasi-parallel bow 


shock. Ion temperature is obtained by calculating the second moment of the velocity 


distribution function and includes the effects of the energetic ions in the foreshock. Panel 
(a) in Figure 2 shows a structure at and upstream of the bow shock consisting of density 
enhancements surrounding a low density region. Examination of panel (a) in Figure 3 
shows the ion temperature in the low density region is over 600 times hotter than the 
pristine solar wind. Note that the ion temperature scale in Figure 3 is set to a maximum 
of 600 for better clarity. This structure looks similar to a simulated HFAs formed at the 


bow shock due to solar wind discontinuities, e.g. Omidi and Sibeck [2007]. Panels (b) 
through (d) in Figures 2 and 3 show the time evolution of this structure that penetrates 
further into the magnetosheath and eventually becomes a part of the highly non-uniform 
and turbulent magnetosheath. In the process the energetic ions within the structure are 
injected into the magnetosheath. 
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To see the signature of this structure and its time evolution as might be observed 
in spacecraft data, Figure 4 shows the ion density, total pressure (normalized to solar 
wind value), velocity (normalized to Va) and temperature, as well as the magnetic field 
(normalized to solar wind value) as observed in time at the location marked by “X” in 
panel (a) of Figure 2. As can be seen, the signature consists of enhancements in density 
and magnetic field (beginning at time -250 Q 1 ) that reach a factor of -3 above the solar 
wind levels. This is followed by large drops in density (minimum value of -15% of solar 
wind density) and field (minimum value of -30% of solar wind magnetic field) in 
association with flow deceleration and deflection and enhancements in ion temperature. 
Note that despite the temperature enhancements, the total pressure in the low density core 


region is below that in the solar wind. Subsequently, the density and magnetic field 


increase above the solar wind levels by a factor of -5 before returning to solar wind 
values. This signature is identical to that of HFAs in general and the SHFAs reported by 
Zhang et al. [2012]. Given the absence of a solar wind discontinuity in the simulation, we 
identify this structure as a SHFA. 


To illustrate the formation of this SHFA, Figure 5 shows the total magnetic field, 
ion temperature and ion velocity in the X direction at two separate times. The top panels 
show a well developed foreshock caviton upstream of the bow shock. The bottom panels 
show that the convection of this caviton by the solar wind into the bow shock transforms 
it into a SHFA. This transformation is associated with further energization of the ions in 
the core of the caviton and the enhancement of the cavity (reduction in magnetic field and 
density) which in turn increases the magnetic field and density in the outer parts. The 
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details of the ion velocity distribution functions within the SHFA and their time evolution 
and their relationship to particle energization process remain to be understood and are 


under investigation. Preliminary results suggest that ion trapping by the cavitons and also 
ion reflection between the bow shock and the cavitons may play an important role in the 


acceleration process. Given the convection of the cavitons towards the bow shock, the 


back and forth motion of ions between the cavitons and the bow shock can result in 


particle acceleration through first and second order Fermi processes 


Examination of the simulation results show that SHFAs form regularly along the 
quasi-parallel bow shock surface as isolated foreshock cavitons, such as that in Figure 5, 
encounter the shock. We also find that at times, multiple cavitons arrive at the bow shock 
near simultaneously and result in the formation of larger and more complex structures. 
An example of this is illustrated in Figure 6 that shows the density zoomed around the 
quasi-parallel shock at 4 different times. Panel (a) in Figure 6 shows the presence of a 
number of SHFA like structures along the bow shock that formed at about the same time 
due to the arrival of multiple foreshock cavitons at the shock. Panels (b) through (d) show 
the time evolution of these SHFAs as they penetrate into the magnetosheath and result in 
large inhomogeneities and turbulence in the quasi-parallel magneto sheath. 

Figure 7 shows the signature of this event in time series data as observed at points “A”, 
“B”, “C” and “D” shown in panel (a) of Figure 6. Density, magnetic field and 
temperature are normalized to solar wind values and flow speed is normalized to the 
Alfven speed in the solar wind. The data looks quite different at each observing point. At 
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point “A”, the data shows signatures associated with 2 SHFAs that are shaded. At point 
“B” two shaded signatures are present that show density and field enhancements and 


depletions, flow deceleration and the presence of energetic ions and look similar to 


SHFAs, however, some differences to SHFAs can also be observed. Similarly, at points 


“C” and “D” signatures similar to SHFAs are present (shaded regions) but clean and full 


signatures of SHFAs are harder to identify. In effect the presence of multiple SHFAs at 
the bow shock and their mutual interactions result in highly nonlinear and complex 
structures whose signatures in spacecraft data would be similarly complex and hard to 
decipher. 


5. SUMMARY AND CONCLUSIONS 

Motivated by the multi-spacecraft THEMIS observations of Spontaneous Hot 
Flow Anomalies at the quasi-parallel bow shock, by Zhang et al. [2012] we have 
examined the structure of a super-critical quasi-parallel bow shock using global hybrid 
simulations. The results show the formation of copious structures at the quasi-parallel 
shock whose time series data resemble those of HFAs and SHFAs. Given the steady 
nature of the solar wind and the absence of a discontinuity in the simulation, these 
structures are identified as SHFAs. The formation of SHFAs in the simulation is tied to 
the convection of foreshock cavitons by the solar wind and their interaction with the bow 
shock. Foreshock cavitons are structures of the order of ~1 Re (Blanco-Cano et al., 2009, 
2011; Kajdic et al., 2010, 2011) consisting of low density and magnetic field core region 
populated with energetic ions and an outer layer with increased density and magnetic 
field strength. Transformation of a caviton to a SHFA is associated with further 
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energization of ions, reductions in density and magnetic field in the core of the cavitons 


and the enhancements of the density and magnetic field in the outer region. The size of 


SHFAs in the Z direction is ~50 ion skin depths which is comparable to that of foreshock 
cavitons and is of the order of 1 Re which is also comparable to the size of HFAs at the 
bow shock. 


Foreshock cavitons have been observed under a wide range of solar wind 
velocities (Mach number) and IMF orientations. During small and intermediate IMF cone 
angles when the foreshock falls upstream of the dayside magnetosphere, foreshock 
cavitons are carried by the solar wind into the bow shock. As a result, we expect the 
formation of SHFAs at the quasi-parallel bow shock over a wide range of solar wind 
conditions. Although the simulation results shown here correspond to Alfven Mach 
number of 12 and IMF cone angle of 10°, examination of other runs with lower Mach 
numbers (down to 6 Va) and cone angles (smaller than 45°) also shows the formation of 
SHFAs at the shock. As such, we believe the formation of SHFAs at the quasi-parallel 
bow shock is a common process and quite significant for ion acceleration and dissipation 
at the super-critical quasi-parallel bow shock. Similarly, the formation and dissipation of 
SHFAs as they interact with the bow shock, is critical for determining the properties of 
the magnetosheath plasma. 

The simulation results also demonstrate that when a number of foreshock cavitons 
arrive and interact with the bow shock near simultaneously, structures larger and more 
complex than SHFAs are formed. These structures are influenced by the interaction of the 
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cavitons with the bow shock but also with each other. As a result, the time series data 


obtained at various points along the bow shock are more complex and varied from point 
to point and exhibit full or partial signatures of multiple SHFAs. Such interactions also 
lead to large inhomogeneities in the magnetosheath. The results presented by Zhang et al. 
[2012] and here demonstrate that ion dissipation processes at the quasi-parallel shock are 
even more complex than previously thought. Future data analysis and simulations are 
needed to shine more light on the impacts of SHFAs on the bow shock, magnetosheath 
and the magnetosphere. Similarly, differences between HFAs and SHFAs and their 


magnetospheric impacts need to be explored further. The fact that the formation of HFAs 


is associated with the presence of solar wind discontinuities while SHFAs form due to the 


interaction of cavitons with the bow shock provide a means of distinguishing between 
HFAs and SHFAs. For example, Zhang et al. [2012] use the absence of a solar wind 
discontinuity associated with an event to identify it as an SHFA. As we leam more about 
SHFAs and how they compare and contrast to HFAs other means of distinguishing 


between the two may become available. 
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FIGURE CAPTIONS 


Figure 1. Panel (a) shows the plasma density normalized to solar wind value and marks 
various parts of the bow shock and the ion foreshock. Panel (b) zooms closer into the 
foreshock and bow shock showing foreshock cavitons. 

Figure 2. Plasma density normali z ed to solar wind value at 4 times (proton gyroperiods 
Q 1 ) demonstrating the interaction of SHFA with the bow shock. 

Figure 3. Ion temperature normalized to solar wind value at 4 times demonstrating 
injection of energetic ions into the magnetosheath by SHFA. 


Figure 4. Time series data showing plasma density, three components of velocity and 
magnetic field and ion temeperature generated at the point marked by “X” in panel (a) of 
Figure 2. Density, total pressure, magnetic field and temperature are normalized to solar 
wind values and velocities are normalized to the Alfven speed in the solar wind. The data 


shows signatures of a SHFA. 


Figure 5. Total magnetic field, ion temperature and velocity in X direction are shown at 
two times demonstrating the transformation of a foreshock caviton into a SHFA. 


Figure 6. Plasma density at 4 times showing the evolution of a number of SHFAs as they 
interact with the bow shock and eventually end up in the magnetosheath. 


Figure 7. Time series data showing the variations of total magnetic field, flow speed 
along X, ion temperature and density at points A, B, C and D marked in panel (a) of 
Figure 6. Density, magnetic field and temperature are normalized to solar wind values 
and flow speed is normalized to the Alfven speed in the solar wind. 
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